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Abstract. Excited bound states are often understood within scattering based 
theories as resulting from the coUision of a particle on a target via a short-range 
potential. We show that the resulting formalism is non-Hermitian and describe the 
Hilbert spaces and metric operator relevant to a correct formulation of such theories. 
The structure and tools employed are the same that have been introduced in current 
works dealing with PT-symmetric and quasi-Hermitian problems. The relevance of the 
non-Hermitian formulation to practical computations is assessed by introducing a non- 
Hermiticity index. We give a numerical example involving scattering by a short-range 
potential in a Coulomb field for which it is seen that even for a small but non-negligible 
non-Hermiticity index the non-Hermitian character of the problem must be taken into 
account. The computation of physical quantities in the relevant Hilbert spaces is also 
discussed. 
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1. Introduction 

The standard formulation of quantum meclianics requires physical observables to be 
mathematically given in terms of Hermitian operators. In recent years theories with 
a non-Hermitian Hamiltonian have been receiving an increasing interest sparked by 
work in the field of PT-symmetric quantum mechanics [T. . PT-symmetric Hamiltonians 
are complex but nevertheless possess a real spectrum. The structure of PT-symmetric 
theories, initially suggested to hinge on the existence of a charge conjugation operator 
[2] has been clarified by showing |2j that the non-Hermitian Hamiltonians could be 
mapped to Hermitian ones, and therefore be fitted within the better-known framework 
of quasi- Hermitian operators Qj. The relevance of the non-Hermitian formulation for 
the description of physical systems is still being debated El E] • 

In the present work we show that the effective Hamiltonians appearing in certain 
theories dealing with bound state scattering by a short-range potential are non- 
Hermitian. In this case the Hamiltonians are real and their non-Hermitian character 
stems from the boundary conditions obeyed by the eigenstates: on the one hand there 
is no physical asymptotic freedom (since the states are bound) and on the other the 
scattering solutions inside the short-range potential region do not exist. However 
contrarily to the situation in PT-symmetric problems, there is in principle an underlying 
Hermitian Hamiltonian, whose solutions are unknown in practice but whose existence 
may provide guiding rules when undertaking practical computations. We will see that 
the unambiguous formulation of bound state scattering sheds some light on issues 
regarding the physical relevance of non-Hermitian formulations of quantum mechanics. 
Let us mention that in the overwhelming majority of applications of the bound state 
scattering formalism to nuclear, atomic or molecular physics non-Hermitian issues have 
been downright ignored; this is unproblematic when non-Hermiticity is small (as is 
generally the case), but we will give an illustration in which ignoring the non-Hermitian 
nature of the scattering Hamiltonian brings in errors that can be directly attributed to 
the (inappropriate) use of the standard inner product. 

We will first introduce bound state scattering theory and show why the scattering 
Hamiltonian is non-Hermitian in the 'physical' Hilbert space Tiph (Sec 2). The quasi- 
Hermitian Hamiltonian will then be described by an expansion in terms of a biorthogonal 
basis, leading naturally to the definition of a new inner product and its associated 
Hilbert space 7i (Sec 3). In line with previous works on quasi- Hemitian operators, we 
will examine the relationship between the two Hilbert spaces Tiph and Ti. in terms of the 
metric operator and further discuss the computation of physical results in TC and Ttph- 
In Sec 4, the formalism will be illustrated by carrying out the numerical calculation of 
an experimentally observable quantity (the autocorrelation function) in the particular 
case of short-range scattering in a Coulomb field. Our concluding remarks will be given 
in Sec 5. 
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2. Scattering description of excited bound states 

Let if^ be the exact Hamiltonian of the 2-particle scattering problem (in the center of 
mass; the physical situation most often considered is that of a light particle colliding on 
a massive compound target). We assume H'^ can be split as 

H' = Ho + V (1) 

where V contains all the short range interactions between the particles. We further 
assume Hq is spherically symmetric (in terms of the relative coordinate) and that short- 
range means that 

V = e{ro-r')Ve{ro-r), (2) 

ie V vanishes outside some small radius r^ {6 is the step function). Therefore Hq 
contains not only the kinetic and internal terms of the non-interacting particles, but 
also any long-range interaction between them. Let E be the total energy; allowing for 
inelastic scattering E is partitioned as 

E = ei + ei (3) 

where and ej are the internal and the kinetic energy respectively (in the case of a 
massive target Si depends on the internal states of the target whereas ej is the collision 
energy of the light particle). The eigenstates of Hq are given by 

\UE)) = \fM)) m)) ; (4) 

fi(ei,r) = (r |/i(ej)) is the eigenfunction of the radial part of Hq whereas the 'target' 
state \i{ei)) includes all the other degrees of freedom, including the non-radial ones of 
the colliding particle (a handy notation given that the angular momenta of the particles 
are usually coupled). The target states are orthogonal, {i\ j) = Sij. Since we are dealing 
with bound states, fi{ei,r) vanishes at and +oo (whenever E is an eigenvalue of Hq). 

The label i defines the scattering channel. In each channel the standing-wave 
solutions are given by the Lippmann-Schwinger equations of scattering theory as 

mE)) = \UE)) + Gq{E)K{E) \UE)) (5) 

where Gq{E) is the principal- value Green's function and K the reaction (scattering) 
operator for standing waves linked to the familiar S matrix by a Cayley transform 
[T^ . The difference here with standard scattering theory is that the bound channels 
are included explicitly [OIIIIIIITT]. The consequences are that (z) Gq{E) has no poles 
- it is modified P relative to the usual resolvent by including a term (solution of the 
homogeneous equation) that has poles at the eigenvalues of H'^ so that overall Gq{E) has 
no poles (but diverges radially) \; (ii) there is no asymptotic freedom: both (r| (pi{E)) 
and (r| ip^lE)) diverge at r — oo for an arbitrary value of E; (Hi) an eigenstate of H*^ 
cannot be given by a single channel solution of the form ^ but requires a superposition 

\r{E)) = J:Z^{E)mE)) (6) 

i 

I As stressed by Fano (Hj who introduced this 'smooth Green's function', for genuine scattering states 
(continuum energies), Go{E) becomes the standard Green's function. 
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where the expansion coefficients Zi{E) are determined by the asymptotic (r oo) 
boundary conditions such that at the eigenvalues (r| il)'^{E)) vanishes at infinity. 

Formally H'^ \'ip'^{E)) = E \iIj'^{E)) is satisfied as well as the usual properties for 
eigenstates of Hermitian operators, such as their orthonormality 

{r{E,)\ r{E2)) = Se,e, (7) 

or the spectral decomposition theorem. However in practice the expansion of Go over 
the eigenstates of Hq is intractable. Instead the radial part of Gq is separated and the 
expansion over the energies reduced to the closed form /j(ei, r^)gi{ei, r>); gi is a solution 
of the radial part of Hq irregular at the origin (for arbitrary bound energies, both / and 
g exponentially diverge in the limit r — > oo). Hence the closed form of the radial Green's 
function only makes sense for r > tq (where V vanishes). This is of course consistent 
with the scattering point of view: when r < r^ we are inside the reaction zone and there 
is no scattering solution, whatever happens within the reaction zone being encoded in 
the phase-shifts. The wavefunction (jH)) outside the reaction zone becomes 



(r mE)) = j:zm 



fi{ei,r) \i) + J29j{^j,r) \j) Kj. 



3 



r > ro 



where Kji are the on-shell elements of the scattering matrix, which are assumed to be 
known. 

It is important to note that the scattering eigenstate (jHl) is the part for r > ro of 
the exact solution \ip^{E)), and not an approximation to it. But within the scattering 
fomulation the 'inner' part of {ip'^lE)) for r < Tq does not exist: all meaningful quantities 
are defined radially on [ro, oo[. Indeed let us write 

\r{E)) = e{r - ro) \^{E)) + 0(ro - r) |V'„(^)) (9) 

and let 

H^Y.E\i^{E)){^{E)\ (10) 

E 

be the restriction of to the outer region r > ro. if is the only operator directly 
known from the solutions of the scattering problem. We have the following properties : 

H = H^ (11) 

(V^(Ei)l V^(E2)) = 6e,E, + ^^iE,E,il - Se,e,) (12) 

H mE)) =^ E \ij{E)) . (13) 

That H is Hermitian relative to the standard product can be seen to follow from its 
definition (fTUj). Eq. (fT^ tells us first that the lip^E)) are normalized to 1 like the \ip'^{E)) 
which might appear surprising in view of © but follows by showing normalization does 
not depend on the inner radial part of the wavefunction (this is done by expressing 
the normalization integral in terms of radial Wronskians, see Sec 5.7 of ,13. ). Eq |T2|l 
also indicates that the scattering eigenstates \ip{E)) are not orthogonal since the scalar 
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product of two disitinct eigenstates is given by /i. This may be shown by rearranging 
eq (jHI) in the form 



where the overall contribution in a given scattering channel i is grouped together. As a 
consequence the radial channel functions Fi{ei,r) must vanish as r — > oo for each i (the 
scattering information is now contained in the F functions and in the new coefficients X 
that both depend on K). Recalling the target states are orthogonal, the scalar product 
(fT^ is seen to depend solely on the radial overlaps between identical channel radial 
functions at different energies, given by 



where W is the Wronskian taken at tq. This equality follows from computing 
(Fj(ei)l |Fj(e2)) — {Fi{e2)\ pi \Fi{ei)) (integrate by parts and recall that the scalar 
product is defined in [ro, oo]). This gives rise to nonzero boundary terms at tq, impliying 
that is not Hermitian on [ro, oo] §. 

Because the \iIj{E)) are not orthogonal, they cannot be eigenstates of the Hermitian 
operator H [eq p3|) ] but are eigenvectors of a non-Hermitian Hamiltonian denoted H. 
From eqs (0) and we see that H is formally given by Hq redefined by restricting 
it radially to the interval [ro, oo] and supplementing it by specific boundary conditions 
on the surface r = ro. It is precisely this fact, combined with the lack of asymptotic 
completeness, that leads to non-Hermiticity. This completes our brief discussion on 
the non-Hermitian character of the bound state scattering problem; we now analyze 
the structure of the non-Hermitian theory and further examine the implications of this 
non-Hermiticity in practical problems. 

3. Quasi-Hermitian operators: metric and Hilbert spaces 

Here we forget about the existence of an underlying exact Hamiltonian and we take 
the practical scattering viewpoint: the phase-shifts are given numbers (obtained from a 
symmetric K matrix) and the physical states are represented by vectors in Tiph, which 
is essentially the Hilbert space of standard quantum mechanics: it is endowed with the 
standard scalar product except that radially the integral is defined on [ro,cxo]. This 
slight modification of the radial integral does not cause any difference since the states of 
interest in scattering phenomena (such as Gaussian states) have negligible probability 
amplitude in the inner zone. In this sense the \ip{E)) belong to Tiph- 

,2 

§ The non-Hcrmitian character of -jpj on bounded intervals with arbitrary boundary conditions is of 
course trivial. In the context of scattering theory, this fact was pointed out in particular by Bloch [T^ 
who introduced a singular surface operator to cancel the boundary terms when defining quantities on 
[0,ro]. However the non-Hermitian character of the scattering eigenstates on [ro,-|-oo] is irrelevant in 
standard scattering theory because the solutions of H'^ and Hq are both (improperly) normalized by 
the same asymptotic condition, hinging on the isometry of the wave operators. 




(14) 




(15) 



£2 - ei 
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Since H is non-Hermitian on Tiph, we have 

{^{E')\H\^{E)) =E{i^{E')\i^{E)) (16) 
mE')\H^ = E' mE') ^E)) . (17) 

We are thus naturally lead to introduce a biorthogonal set { 'ip{E)^ , \%1){E))} 17J, where 
we denote by iIj{E)) the eigenstates of . The following properties are satisfied: 

(18) 
(19) 

(20) 



H\ij{E))=E\i;{E)) 



HE) 



E 



[i;{E) ^E')) = 6ee' 
from which it follows that we can write the following expansions: 

H = J2E im)) {i'iE)\ H^ = Y.E \i'{E)) (^(E)| 

E E 

H and are further linked by 

where ^ is a Hermitian operator given by 

G = j:\he)) {^{E)\ 



E 



(21) 

(22) 

(23) 
(24) 



We will take for granted the completeness of the biorthogonal basis, although it is by no 
means obvious. In particular the difficulties that arise when Tiph is of infinite dimensions 
have been pointed out recently [T31 Hn] . Completeness of the biorthogonal basis implies 
that the 'canonical metric basis' (in the sense of [Sj), consisting of the eigenvectors of the 
metric operator, is also complete. From there we deduce that an arbitrary state of TCph 
can in principle be expanded in terms of the \ijj{E)), ie the eigenstates of H span the 
entire Hilbert space of admissible physical states even if they do not form an orthogonal 
basis in Hph- 

The relations p6|) - P^ have become familiar lately in the context of PT-symmetric 
quantum mechanics and more largely in works dealing with quasi-Hermitian operators 
(see in particular [5]). Eq (j22I) is the defining relation of quasi- Hermiticity jTHI provided 
g is invertible {g~^ then being its inverse, since by ()20j) gg~^ is a representation of the 
unit operator in H) and positive-definite. We will not attempt to prove these properties 
here. We note however that if the \iIj{E)) (and hence the il){E)^) form a basis of 7ip/i, as 
we have assumed to be the case, then g has no null eigenvalue and is thus invertible. It 
is of course a working hypothesis in scattering theory that any meaningful physical state 
can be expanded in terms of the \il){E)) (but this may not be true mathematically for a 
given arbitrary vector). The positive-definiteness of g follows heuristically by remarking 
that in the 'mixed' representation 



^(^0 g\^{E)) 



(25) 
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simply becomes (|12p . so that Q > I + where I is the identity matrix, M is the 
special matrix with elements Mij = 1 — 5ij and /i a small (|yu| <^ 1) real number. The 
positive-definiteness of / + /iM ensures that Q is positive definite too. The positive- 
definitiness of Q is important to define a positive norm in 7f jH ^j. Since from eq 

(121 

^{E)) = \i,{E)) , (26) 
the inner product is defined through 

(^(Ei), V^(E2))g = {i^{E,)\G-^ |^(E2)) = {^{E,) |^(E2)) = 6e,e,. (27) 

Q is thus seen to be (the positive definite) metric. By eq (j^^ it is immediate to verify 
that H is Hermitian relative to this new inner product. 

Let Ti be the Hilbert space endowed with the inner product defined by 
Calculations are simple to perform in Ti because the new scalar product reestablishes 
orthogonality. Indeed let |0i) = Y,Cii{E)\ip{E)) and |02) = Y,Ci2{E)\ilj{E)) be two 
vectors in Ti. Then it follows from eq ipTj) that 

(01, 01)^ = (0i|0i) = ^|ai(E)|2 = l (28) 

E 

(01,02)g = (01 102) = E«K^)«2(i?) (29) 

E 

with the obvious notation 

^Q-'\ct^i) = Y.^i{E)g'^\i,{E)). (30) 

E 

We further see that quantities involving the expansions of the non-Hermitian 
Hamiltonian, such as the time evolution operator, cannot be directly determined in 
H-ph-, precisely because of the non-Hermiticity of H in Tiph- But in Ti the evolution 
operator is given by 

U{t) = Y.e-''''\m)){i'{E)\. (31) 

E 

Hence for example if we take an initial state as |0(t = 0)) = |0i), the state evolves 
according to 

\m = E e"^""* (V^(E), 00, = E e-''''a,{E) |^(E)) , (32) 

E E 

operating in effect in Ti as we would in Tiph with a Hermitian operator. 

However, in scattering problems, the physical states are known in Tiph-, not in Ti. 
Let ICi) and IC2) be two vectors in Tiph and assume they can be expanded over the 
\%Ij{E)) as ICi) = Y.ciii.E) \iij{E)). They are normalized relative to the standard scalar 
product, 

(0 10 = 1 = E <{E)a,{E') {ij{E) mE')) ; (33) 

EE' 

since the basis is nonorthogonal in Hph, J2e \(^i{E)f 7^ 1- On the other hand operators 
involving the Hamiltonian, such as the evolution operator (p?T|l are known on 7i but not 
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on Tiph- The transformation between the two Hilbert spaces can be done both ways, 
H — > Hph for the operators or Hph — ^ Ti. for the states. Indeed if an operator A is 
Hermitian in Ti then 

A = g-^/^Ag^/^ (34) 

is Hermitian in Tiph- This folfows directly from the general version of eq ()22j) . 

A = gA+g-\ (35) 

This transformation defines a linear map ^ that leaves the inner product invariant: 

(01,02)g = (01 102) = {<Pl\S-' 102) = (Cll C2) (36) 

where we have defined 

10) = g-'/' 10.) . (37) 

Therefore \Q) and represent the same physical state but relative to different Hilbert 
spaces: IQ) in Hph and in H. Of course as vectors we may as well have for instance 
\(f)i) G Hph but then does not describe the same physical state as it does in H. 
It is interesting to note that the functions \iIj'^{E)) defined on Hph (with r G [0,+oo]) 
represent the exact eigenstates of the underlying Hamiltonian. But the lip^E)) envisaged 
as the restriction for r > ro of the \ip'^{E)) do not represent the eigenstates in Hph (now 
with r G [ro, +00]) but in H, that is on the Hilbert space in which the Hamiltonian H 
is Hermitian, despite the fact that (r \iIj{E)) and (r |?/''^(i?)) are identical for r > tq. 

Finally, we briefly describe how to undertake practical calculations. Recalling that 
scattering solutions as well as the Hamiltonian H are defined in H, and comparing eqs 
()34|) and (jHTj) . it appears that it is computationally simpler to transform the physical 
states from Hph to H rather than transform the operators to Hph- Nevertheless in both 
cases it is necessary to determine the metric g. In general (as in the illustration given 
below) ^ is a matrix of infinite rank=. g is therefore truncated around the energy 
interval of interest. The matrix elements are determined in the 'mixed' representation 
given by eq which simply amounts to determine the overlaps 

gEE' = mE)\^iE')) (38) 

where E and E' span the (truncated) finite interval. The resulting matrix g~^ is 
numerically inverted, allowing to determine the second set of the biorthogonal basis 
by eq (j^ . g can also be diagonalized, retrieving in a single step g~^, g~^^'^ and g^^'^; 
we then compute the operators in Hph or the representation of the physical states in H 
by inverting eq (jSZ|)- The degree of non-Hermiticity is assessed through the metric in the 
mixed representation ()38|): if the Hamiltonian is Hermitian relative to the standard inner 
product, g becomes the identity matrix. As non-Hermiticity becomes important, the off- 
diagonal elements of the metric increase. To assess the degree of non-Hermiticity we 
introduce a non-Hermiticity index k that we define somewhat arbitarily by the average 
of the N largest absolute values of ^ — J (ie the N largest off-diagonal terms of the 
metric) where N is the dimension of the chunk of g under study, k is thus a local 
spectral measure of non-Hermiticity. 
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Figure 1. Each dot represents one of the 200 {N/2) largest off-diagonal elements 
^BB'I of the metric. All the diagonal elements are given by Qee = 1- 



4. Illustration 

To illustrate the formalism given above we will take an example in the context of the 
bound states formed by the scattering of an electron on a positively charged target. This 
situation is widely employed in atomic physics to study the highly excited ('Rydberg') 
states of atoms with a single excited electron. More specifically we will compute the 
autocorrelation function 



in two ways: by taking into account the non-Hermitian character of the Hamiltonian on 
the one hand, and by downright ignoring Hermiticity related issues on the other hand. 
|C(t)| is in principle an experimentally observable quantity. If the no n- Hermiticity index 
K is negligible, the two methods of calculation will give nearly identical result (for typical 
atoms K turns out to be very small, although non-Hermitian issues have always been 
ignored from first principles). 

The long-range Hamiltonian Hq in eq (0) contains the radial Hamiltonian of the 
colliding electron in a centrifugal Coulomb potential as well as the free Hamiltonian of 
the target (an atomic ion). /j(ej, r) in eq Q is therefore a Coulomb function regular at 
the origin (it is also regular at +00 only if e belongs to the spectrum of the radial part 
of Hq, ie when e = — l/2n^, n G A^). The radial channel functions Fi{ei, r) appearing in 
eq (HH), solutions of the radial part of the redefined Hq for r > tq, are given by a linear 
combination of Coulomb functions regular and irregular at the origin, the combination 
ensuring that Fj(ej,r) converges at 00 ||. For the scattering matrix K{E) we take a 
6x6 matrix with a strong energy dependence. We also set the 6 values of Ei to model 
the internal energies of the target (we take ei = for the ground state and 5 different 
values for the excited states of the target). The bound state energies and coefficients 

II Note that Fi{ei,r) mathematically diverges as r ^ 0, which is of course irrelevant to the scattering 
problem defined on [ro, +00] 



C(t) = (C(t = o)|C(t)) 



(39) 
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Figure 2. The autocorrelation function C(t) is given as a function of time (in units 
of the period of the classical orbit of the electron having the mean energy of the initial 
state). The inset shows the short-time situation. Dashed (colour online: red) curve: 
the autocorrelation function is computed by ignoring the non-Hermitian character of 
the problem, following eq 143|l . Solid line: Result computed by taking non-Hermiticity 
into account (eq (0HJ)). 




are obtained by applying the boundary condition (r \%Ij{E)) — ^> as r — > oo, yielding the 
system [TH] 

[K{E) + R{E)] Z{E) = (40) 

where R{E) is a diagonal matrix with elements R{E)ii = tan7r(— 2(£' — This 
system is solved numerically for E and then the nontrivial solutions Zi{E) are obtained. 
We compute about N = 400 eigenstates. The radial overlaps (fT3j) are determined 
analytically, and from there we compute the metric elements Qee'- For the overall 
chunk, the non-Hermiticity index is calculated as k = 0.07. The ordered distribution of 
the N largest off-diagonal elements of the metric is shown in Fig. 1. 

We now choose an initial state \C(t = 0)), that we take to be a Gaussian localized 
radially very far from the target, at the outer turning point of the radial potential for an 
excited electron (with a mean energy n = 55), with the target being in its ground state. 
Initially \((t = 0)) is defined on an orthogonal basis of Hph but we assume (and verify 
numerically) that this state can approximately be expanded on our chunk of computed 
eigenstates of H as 

\at = o)) = j:a{E)mE)). (4i) 

where the a{E) are determined by projections. At this point we proceed along the 
two different lines mentioned above. In the first method of calculation we employ 
the machinery of standard (Hermitian) quantum mechanics, ignoring non-Hermiticity 
issues. This may appear absurd in view of the preceding discussion, but this is the way 
computations are undertaken in applied problems ^. Moreover this will allow us to 

% It is true that in typical atomic problems, k is significantly smaller (below 10^'^) than in the example 
given here, so that the computed results would only be marginally affected by taking into account the 
non-Hermitian character of the Hamiltonian. 
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assess the relevance of the formahsm given above in practical calculations - as we will 
see by comparing the first method to the second one, where the formalism developed in 
Sec 3 will be employed. 

In the first method the expansions 

J2 ItPiE)) {ipiE)\ [1 or E or exp{-iEt/n)] (42) 

E 

are taken as representations of the unit operator, the Hamiltonian or the evolution 
operator respectively. The coefficients a{E) of eq (PT|) are thus given by the projection 
of this 'unit' operator as {ip{E) \ ({t = 0)), and the autocorrelation function (|39|l follows 
by employing this 'evolution' operator, 

C(t) = 5:e-^*/^|WE)|C(t = 0))r. (43) 

E 

The result is shown in Fig. 2 by the dashed line; in particular the inset shows the short- 
time evolution, and it may be noticed that at t = we do not have C(t = 0) = 1, ie 
|C(t = 0)) is not normalized after the application of the 'unit' operator which as 
we know is not the correct unit operator on Tiph- Neither is the 'evolution' defined by 
eq (jl^ unitary: {((t) \ ((t)) computed with eq (jl21) shows strong oscillations, displayed 
in Fig. 3. 

The correct method to compute C{t) involves first mapping \({t = 0)) to 7i, yielding 
^^/^ \((t = 0)) (cf eq dnZI)), then apply the unitary evolution operator in H given by eq 
()31|) and finally compute the result with the inner product (jHEI) in 'H. If we follow the 



notation (jHTj) and put 

m = 0)) = g'/' \at = 0)) (44) 

we get the following equivalent expressions for the autocorrelation function: 

C(t)= (0(t = O),t7(t)0(t = O))^ (45) 

= (0(t = O)|f7(t)|</)(t = O)) (46) 

= {<p{t = o)\g-'u{t)m = o)) (47) 

= {at = o)\g-'^'uit)g'/'\at = o)). (48) 



Eqs fl45j) and ()4(j|l give the autocorrelation function as computed entirely in Ti. whereas 
eq fliH|) is the same expression in Hph- g^^^'^U{t)g^^'^ appears as the (correct and 
unitary) evolution operator in Hph resulting from the mapping given by eq (jUj). The 
computed result is shown by the solid line in Fig. 2, which of course obeys C(t = 0) = 1 
(normalization at other times follows from unitarity). 

The most salient feature arising from the comparison of the two curves in Fig. 2 
concerns the different profiles of the autocorrelation functions. This implies that it will 
not be possible to recover the correct result from the first method result ()43j) by 
simply renormalizing the latter in Tiph (as is sometimes done in practical scattering 
problems). Conversely it would not make much sense to assume that the initial physical 
state pT|) is known in H, so that one would not need to determine mapping H Hph- 
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Figure 3. The 'norm' {C,{t)\ C{t)) is computed as function of time (in units of the 
period of the mean Kepler orbit). The norm is not conserved because the 'evohition' 
operator 142|l is Hermitian but not unitary in Hph- 



Such an exception happens in the specific but nevertheless important cases in which one 
is only interested in transitions involving given eigenstates of H. 

5. Conclusion 

We have seen that the widely employed formalism of bound state scattering theory 
should be properly understood within the framework of non-Hermitian quantum 
mechanics. Although in typical cases the non-Hermiticity index n is small so that in 
practice non-Hermitian issues can be ignored, we have given an illustration for which the 
calculations of experimentally observable quantities require the proper non-Hermitian 
formulation. The latter has essentially the same structure and tools as the PT-symmetric 
systems (reformulated in the quasi-Hermitian framework) that are currently being 
extensively investigated. However in the present case, the physical meaning of non- 
Hermiticity is more transparent than in the case of PT-symmetric quantum-mechanics. 
In particular, we have seen that changing the radial interval minutely from [0, +oo] 
in the underlying exact problem to [ro, -|-cxd] in the scattering problem leads to an 
entire redefinition of the Hilbert spaces relevant for quantum mechanics. Indeed, by 
this change the Hamiltonian H becomes quasi-Hermitian on Tiph- One can then either 
redefine the inner product, constructing a new Hilbert space H, or map the states 
and operators to the physical Hilbert space Hpt- We have seen that computations are 
simpler to undertake in Ti than in Tiph, but except in the specific cases involving the 
sole eigenstates of the non-Hermitian Hamiltonian, this simplicity is only apparent: as 
arbitrary physical states are known in Hph, the mapping between the two Hilbert spaces 
must be explicitly determined anyway, involving the computation of the metric. In the 
example given in this work the metric was constructed from the numerical calculation 
of the exact eigenstates of if in a restricted energy interval of interest. 

The interesting insight gained by the existence of an underlying exact Hamiltonian 
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H'^ is that for scattering states, the physical Hilbert space Hph is essentially the same as 
the Hilbert space of the exact problem. But the expansion in Tiph of a physical state in 
terms of the eigenstates of the exact Hamiltonian differs from the expansion in terms of 
the eigenstates ipiE) of the non-Hermitian Hamiltonian (although the physical results 
- eigenvalues, probability amplitudes - will be identical). Actually the expansion in 
terms of the eigenstates of the exact Hamiltonian in Tiph is identical to the expansion in 
terms of the eigenstates of H in Ti. These remarks suggest that as far as the scattering 
eigenstates are concerned, H is more physical than Hph, where these eigenstates become 
\^(E)). From a more general standpoint it appears that quantum mechanics 
requires above all a Hilbert space Ti. on which the operators are self-adjoint relative to 
a given inner product, whatever this inner product may be. In this work the Hilbert 
space Hph defined with the standard inner product (the inner product) only entered 
the problem because in bound state scattering arbitrary physical states and operators 
are already known in this space. In general however it is possible to envisage the case 
in which the standard inner product would not play a special role, although such a 
situation will probably lead to intricate interpretational issues regarding the physical 
significance of computed quantities. 
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